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Summary 


•  In  Primary  Auditory  Cortex  (AI)  of  Ferrets,  we  have 
previously  characterized  cells’  responses  to  dynamic 
broad-band  sounds.  We  found  best  responses  to 
temporal  modulations  from  4  to  16  Hz,  and  spectral 
modulations  from  0.4  to  1 .6  cycles/octave  in  the 
stimulus’s  spectro-temporal  envelope. 

•  The  Spectro-Temporal  Response  Field  (STRF) 
explains  the  linear  component  of  the  response  to  the 
spectro-temporal  envelope  of  a  broad-band  sound. 

•  The  STRF  is  often  a  good  predictor  of  the  response  to 
an  arbitrary  sound.  However,  previous  measurements 
of  the  STRF  using  sinusoidal  spectro-temporal 
envelopes  were  hampered  by  the  time  required  to 
accumulate  data  from  a  cell. 

•  We  use  sums  of  spectro-temporal  sinusoids  as  stimuli: 
these 

-  reduce  recording  time 

-  confirm  quadrant  spectro-temporal  separability 

-  can  be  used  to  explore  non-linearities 
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Ripples,  or  Auditory  Gratings 

Ripples  are  auditory  gratings  whose  spectral  envelope  is  a  sinusoid  along 
the  log(frequency)  axis.  At  any  time  t  and  any  frequency  x,  the  amplitude 
is  given  by: 


5(^)=  sin[2nwt  +  2nClx  +  (|)]  w  -  ripple  velocity 

x  =  log 2\f/f0]  Cl  =  ripple  density 
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The  Ripple  Domain  is  the  Fourier 
space  of  the  spectrograms.  We 
probe  a  cell  at  different  velocities 
w  and  different  densities  Cl ,  and 
quantify  the  response  for  up  and 
down-moving  sounds. 

Any  ripple  in  the  lower  half-plane 
is  equivalent  to  a  ripple  in  the 
upper-half  plane. 
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Spectro-Temporal  Noise 

In  order  to  speed  up  the  characterization  of  a  cell’s  response,  we  have  used 
multiple  combinations  of  ripples,  of  all  velocities  w  and  densities  £ 1 ,  with 
random  phases.  Different  combinations  have  different  choices  of  individual 
ripple  phases.  The  range  of  frequencies  and/or  velocities  is  adjusted  to  match 
the  range  of  interest  to  the  cells  being  studied.  We  generally  use  -24  Hz  to  24 
Hz  for  cortex,  and  -400  Hz  to  400  Hz  for  the  Inferior  Colliculus. 

This  is  a  Spectro-Temporal  generalization  of  white  noise. 


,S'noise«rv)  =Yjj  Tjk  sin[2jw/  +  InQpc  + 


Examples: 

Cortical  Stimulus 


Mid-Brain  Stimulus 


C/^A/2  1 

Center  for  Auditory 
and  Acoustic  Research 


|_  _ M M M MMM MM M M M M I^H 

-400  0 


w 


400 


Institute  for  Systems  Research 
University  of  Ma  tyla  nd 


STRF  from  Cross-Correlation  with  Noise 


Stimulus:  S(t,x) 

Response:  R(t )  =  ^  5(f  -  tk ) 

*  i  rT 

Cross-Correlation:  C(t,  x)  =  x)R(f-x)dt 

|  Jo 


(=  Spike-Triggered  Average) 


C(x,  x)  contains  cross  terms 

The  cross  terms  have  random  phase  and  can  be  attenuated  by 
averaging  over  multiple,  random-phase  stimuli. 
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Temporally  Orthogonal  Ripple  Combinations 

(TORCs) 

•  Stimuli  are  composed  only  of  ripples  with  different  ripple  velocities. 

•  Each  stimulus  contains  ripples  which  cover  the  same  range  of  ripple 
velocities,  but  at  different  ripple  frequencies. 

•  Multiple  stimuli  are  still  needed  to  present  a  complete  set  of  ripples. 
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STRF  from  Cross-Correlation  with  TORCs 


•  TORCs  are  better  suited  for  temporal  cross-correlation  because  there 
are  no  cross  terms. 

•  The  resulting  estimates  are  robust,  use  short-duration  stimuli,  and  are 
quickly  computed. 


strf™rc 
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Simulation 
STRFest(x,  x) 
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TORC  Cross-Correlation  vs.  Single-Ripple 


Single  Ripple  Cross-Correlation  Single  Ripple  Cross-Correlation 


Center  for  Auditory 
and  Acoustic  Research 


Institute  for  Systems  Research 
University  of  Ma  tyla  nd 


/(kHz) 


Predicting  Responses  from  STRF 

The  response  to  an  arbitrary  sound  is  given  by  the  convolution  of  the 
STRF  with  the  stimulus  envelope,  plus  a  constant. 


=  X  L  *)  *'  S(f,  x)}  +  E{/f(f)} 
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Stimuli  Used  for  Predictions 
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Predictions  and  Responses  I 


STRFest(f,  x) 
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Predictions  and  Responses  II 
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Non-Linearity 

•  The  value  of  the  STRF  at  each  point  is  the  slope  of  a  linear  rate-level  function: 

R  ( t )  =  STRF(x,  x)  ■  S(t-% ,  x)  .  The  instantaneous  rate  is  given  by  a  linear  operation  between  the 
stimulus  envelope  and  the  STRF. 

•  Polynomial  rate-level  curves  can  be  measured  at  every  (x,  x),  and  improve  the  description. 

•  The  coefficients  of  a  power  series  expansion  of  the  curves  are  given  by  the  diagonals  of  the 
Volterra  kernels. 

•  Using  cubic  polynomials,  and  inverse-repeat  stimuli,  we  have  shown  that  either  the  non- 
linearities  are  absent,  or  they  are  restricted  to  second  order. 

•  Subtraction  of  the  response  to  the  inverted  envelope  from  the  response  to  the  non-inverted 
envelope  gives  a  polynomial  fit  dominated  by  the  linear  term.  This  would  be  expected,  for 
example,  from  a  rectifying  non-linearity. 


-  -  -  -  Measured  rate-level  function  at  one  x  and  v 

-  3rd  Order  Regression  Curve 

-  3rd  Order  Regression  with  Inverse-Repeat 

-  Linear  Transfer  Characteristic  given  by  STRF  Estimate 

-  Mean  Spike  Rate 
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Spectro-Temporal  Rate-Level  Functions 


Rate-level  functions  change  with  x  and  v. 


<D 


Normalized  Stimulus  Level 
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Non-linearity  Theory 


The  STRF  is  the  linear  time-invariant  filter  governing  the  transformation 
from  stimulus  to  response.  Thus  it  can  be  identified  with  the  first  kernel 
of  a  Volterra  series  expansion  of  the  system. 


Third-order  Volterra  series  expansion  of  an  input  S(f,  x): 

R(t)  =  v0  + 1 dz^ dx  •  vl(z,x)S(t  -  z,x) 

+| dz1  J  Jt2| dxl  J dx2  •  v2(zl,z2,xl,x2)S(t  -  zvx2)S(t  -  z2,x2) 

+J dzx  J dz2^dz3  J dxx  J dx2  J dx3  •  v3(zvz2,z3,xvx2,x3)  • 

S(t-zl,xl)S(t-z2,x2)S(t-z3,x3)  +  ••• 


Form  of  the  Regression  Function: 
(Third-order  Approximation) 

g(s;z,x)  =  F{R(f)  I  S(t  -z,x)  =  s } 

~  a0(z, x)  +  ax(z,x)  •  s  +  a2(z,x)  •  s 2  +  a3(z, x)  •  s3 


•  The  regression  functions  describe  the  non-linearities  within  each 
channel,  but  not  interactions  between  channels. 

The  true  STRF  is  just  the  linear  part... 

STRF(z,  x) = ax  (t,  v) 


Use  of  the  Regression  Function 
to  Predict  the  Response: 

R(t)  =  | dz | dx  •  a0(z,x)  +  J dz J dx  •  ax(z,x)S(t  -  z,x) 

+|  dr  J  dx  •  a2  (t,  v)S2  (t  -  z,  x)  +  J  dr  J  dx  •  a3  (z,  x)S 3  ( t  -  z,  v) 
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